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Abstract 

We present the result of a direct search for the decay Ks — > e + e~ , obtained with a 
sample of e + e~ — ► <p — > K$Kl events produced at DA$NE, the Frascati ^-factory, 
for an integrated luminosity of 1.9 fb . The search has been performed using a 
pure Ks beam tagged by the simultaneous detection of a Kl interaction in the 
calorimeter. Background rejection has been optimized by using both kinematic and 
particle identification cuts. We find BR(Ks — ► e + e~) < 9 x 10~ 9 at 90% CL, which 
improves by an order of magnitude on the previous best limit. 

Key words: e + e~ collisions, DA<I>NE, KLOE, rare Ks decays, CP, xPT 
PACS: 13.20.Eb 



1 Introduction 

The Kl, Ks decays into leptons pairs (e + e~, are due to AS = 1 

flavour-changing neutral-currents (FCNC) transitions. The decay amplitudes 
receive contributions both from long distance (LD) effects, dominated by the 
27 intermediate state shown in Fig. 1, and from short-distance (SD) effects, 
due to box and penguin diagrams via W,Z exchange. The SD contribution can 
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Fig. 1. Long distance contribution to Ks — > £ + £ process, mediated by two-photon 
exchange. 

be rather precisely evaluated in the Standard Model (SM) so that a theoretical 
evaluation of the LD contribution would allow a comparison of experimental 
results against predictions that may represent a significant test of the SM. For 
the Kl decay the evaluation of the LD contribution needs a model for the 
Kl — » 7*7* form factor, while for the Ks it can be determined at lowest order 
in the chiral perturbation theory. In this case one obtains [1]: 

T(K S - n+n-)/T(K s -> 77) *l 2 x 10" 6 , 
T(K S -> e + e-)/r(A^ -> 77) ~ 8 x 1CT 9 . 
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Using the present average [2] for BR(J^5 — > 77), we evaluate BR(Ks — > 
e + e~) ~ 2x 10~ 14 . A value significantly higher would point to new physics. The 
best experimental limit for BK(Ks — > e + e~) has been obtained by CPLEAR [3], 
and it is equal to 1.4 x 10 -7 , at 90% CL. Here we present a new search for 
this decay, which improves on the previous limit by more than an order of 
magnitude. 



2 Experimental setup 



The data were collected with the KLOE detector at DA$NE, the Frascati 0- 
factory. DA$NE is an e + e~ collider that operates at a center-of-mass energy 
of ~ 1020 MeV, the mass of the <fi meson. Positron and electron beams of equal 
energy collide at an angle of (tt — 0.025) rad, producing mesons with a small 
momentum in the horizontal plane: ~ 13MeV/c. mesons decay ~ 34% 
of the time into nearly collinear K°K° pairs. Because J PC (<f)) = 1 , the 
kaon pair is in an antisymmetric state, so that the final state is always K$Kl. 
The contamination from K^Kl and K$Ks final states is negligible. Therefore, 
the detection of a Kl signals the presence of a K$ of known momentum and 
direction, independently of its decay mode. This technique is called K$ tagging. 
The sample analyzed corresponds to an integrated luminosity of ~ 1.9 fb -1 , 
yielding ~ 2 billion KqK^ pairs. 
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Fig. 2. Vertical cross section of the KLOE detector. 



The KLOE detector (Fig. 2) consists of a large cylindrical drift chamber (DC), 
surrounded by a lead/scintillating-fiber sampling calorimeter (EMC). A super- 
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conducting coil surrounding the calorimeter provides a 0.52 T magnetic field. 
The drift chamber [4] is 4 m in diameter and 3.3 m long. The chamber shell 
is made of carbon-fiber/epoxy composite, and the gas used is a 90% helium, 
10% isobutane mixture. These features maximize transparency to photons and 
reduce Kl^Ks regeneration and multiple scattering. The momentum resolu- 
tion is a(p±)/p± = 0.4%, and the Kg — > ir + 7i" invariant mass is reconstructed 
with a resolution of ~ 1 MeV. 

The calorimeter [5] is divided into a barrel and two endcaps, covering ~ 98% 
of the solid angle. The modules are read out at both ends by photomultiplier 
tubes. The arrival times of particles and the three-dimensional positions of the 
energy deposits are determined from the signals at the two ends. The read- 
out granularity is ~ 4.4 x 4.4 cm 2 ; fired "cells" close in space and time are 
grouped into a "calorimeter cluster". For each cluster, the energy E c \ is the 
sum of the cell energies, and the time t c \ and the position x c \ are calculated as 
energy-weighted averages over the fired cells. The energy and time resolutions 
are a E /E = 5.7%/ v / J E(GeV) and a t = 57ps/ v / / £(GeV) © 100 ps, respectively. 

The trigger [6] used for this analysis requires two local energy deposits above 
a threshold of 50 MeV in the barrel and 150 MeV in the endcaps. The trigger 
has a large time spread with respect to the beam crossing time. However, it 
is synchronized with the machine RF divided by 4, T sync ~ 10.8 ns, with an 
accuracy of 50 ps. An estimate of the event production time (T ) is determined 
offline. 

The response of the detector to the decays of interest and the various back- 
grounds were studied by using the KLOE Monte Carlo (MC) simulation pro- 
gram [7]. Changes in the machine operation and background conditions are 
simulated on a run-by-run basis. The beam energies and the crossing angle are 
obtained from the analysis of Bhabha scattering events with polar angle 
above 45 degrees. The average value of the center-of-mass energy is evaluated 
with a precision of 30keV for each 100nb _1 of integrated luminosity. 

To study the background rejection, a MC sample of </> decays to all possible fi- 
nal states has been used, equivalent to an integrated luminosity of ~ 2.1 fb -1 . 
A MC sample of 45 000 K$ — > e + e~ events has been also produced, corre- 
sponding to a BR of 1.6 x 10 -4 . This sample is used to measure the selection 
efficiency, and includes radiative corrections. Two processes are expected to 
contribute to photon emission: the inner bremsstrahlung photon emission, 
Ks — > e + e~ + a two-photon decay with one virtual photon conversion, 
Ks — > 77* — > ■ye + e~. The first process is simulated using the PHOTOS [8] 
generator. The events due to the second process are rejected by the kinematic 
cuts used in the analysis. 
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3 Data analysis 



The identification of i^-interaction in the EMC is used to tag the presence of 
Ks mesons. The mean decay lengths of K$ and Kl are As ~ 0.6 cm and A^ ~ 
350 cm, respectively. About 50% of K^s therefore reach the calorimeter before 
decaying. The Kl interaction in the calorimeter barrel (-K'crash) is identified by 
requiring a cluster of energy above 125 MeV, not associated with any track, and 
with a time corresponding to the Kl velocity in the <ft rest frame, (3* ~ 0.216. 
Requiring 0.17 < (3* < 0.28 we selected ~ 650 million i^-tagged events 
(i^ cras h events in the following), which are used as a starting sample for the 
Ks — > e + e~ search. 

As a first step of the signal search, we select events with tracks of opposite 
charge having point of closest approach to the origin within a cylinder 4 cm in 
radius and 10 cm in length along the beam line. The two tracks are required 
to form a vertex with position in the transverse plane p < 4 cm. Moreover, the 
track momenta (p) and polar angles (9) must satisfy the cuts: 120MeV/c < 
p < 350MeV/c and 30° < 9 < 150°. The tracks must also reach the EMC 
without spiralling, and have an associated cluster with energy E c i > 50 MeV 
and position in the transverse plane p c \ > 60 cm. These requirements ensure a 
redundant determination of the event-To and allow us to evaluate the time of 
flight (TOF) for each particle. 

The two-track invariant mass evaluated in the hypothesis of electron mass, 
M ee , is used to reject the dominant background due to Ks — > We 
require M ee > 420 MeV, thus removing most of Ks — > tt + 7t~ events which 
peak at M ee ~ 409 MeV, with a resolution of ~ 1 MeV. In order to reject 
tracks with a larger uncertainty on the fit parameters, we also require the 
propagated error on the invariant mass, AM ee , to be less than 2.5 MeV. In 
Fig. 3, the M ee distribution is shown for both MC signal and background. 
The background is due to the following sources: residual Ks — > tt + tc~ events, 
populating the low M ee region, and cb — > 7r + 7r~7r° events, spreading over the 
whole spectrum. A Ks —> 7r + 7r~ event can have such a high value of M ee either 
because one track is badly reconstructed {ttti component in the following) or 
because one pion decays to a muon before entering the DC and the vertex 
is reconstructed from a pion and a muon track (irp component hereafter). A 
(p — > 7T + 7r _ 7r° event (3ir component in the following) can satisfy the Ks tagging 
criteria from the presence of a machine background cluster (fake i^ cra sh)- At 
this stage of the analysis, we are left with ~5x 10 5 events. The efficiency for 
signal selection, given the Ks tagging, is ~ 0.54, as evaluated using MC. 

The absolute background level is not taken directly from MC , but is obtained 
by normalization of data in the region of signal sidebands. The reliability of 
MC prediction is checked comparing with data after each step of the analysis. 
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Fig. 3. M ee MC spectra for signal (open histogram) and the main background 
sources (gray histograms), as explained in the text; the signal corresponds to a 
BR of 1.6 x 10" 4 . 

For this purpose, the M ee interval is divided into a signal region, around the 
kaon mass peak, and two sidebands: 

420 < M ee < 460 MeV region 1 (left sideband), 

460 < M ee < 530 MeV region 2 (signal), (2) 

530 < M ee < 700 MeV region 3 (right sideband). 



7T7T and 7r/x background sources largely dominate on 3tt component in region 1, 
the opposite occurring in region 3. A scale factor for the 3ir component, / 37r , 
is therefore directly evaluated in region 3 as the ratio of the number of events 
in the data sample and the number of MC 3n events. We obtain: 



which has to be compared with a data/MC luminosity scale factor of ~ 0.86. 
The observed discrepancy is well understood, being due to the fact that MC 
underestimates the rate of fake -K" cras h from machine background. After nor- 
malization, the M ee shape is well reproduced by MC 3n sample, as shown in 
the left panel of Fig. 4. To obtain the scale factors f nn and f vll for the 7T7t and 
7r/i components, we fit the M ee distribution in region 1 to a linear combination 
of the MC background spectra with the 37r component fixed as in Eq. 3. The 
MC distribution after fit is compared to data in the right panel of Fig. 4; we 
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obtain: 

= 0.861 ±0.005, 

= 1.249 ± 0.008, (4) 
p = -0.77, 

where p is the correlation factor and the errors quoted are statistical only. The 
scale factors f nfJi and f nn have to be compared with the expected data/MC 
ratio of ~ 0.73. A sizable deviation is observed for tttt events, which is expected 
since MC tends to underestimate the rate of events in the very far tails of the 
tracking resolution. After normalization, the number of background events, 
, is estimated as: 

N™ = U x + U x K° + f** >< <° ( 5 ) 



> 1000 



- 800 
a 

> 

W 600 



400 - 



200 



MC 3ji 
Data 



++ + + 



++ 



MC Tot Bkg 

■ MC TC7C 

■ MC tc(J, 
Data 



540 560 580 600 620 640 660 680 700 

M„ (MeV) 




M (MeV) 



Fig. 4. Data-MC comparison for M ee spectra in region 3 (left) and region 1 (right), 
after normalization; data are represented by the black points, MC background com- 
ponents by gray histograms. 



Kinematics and topology can be further exploited to improve background 
rejection. For most of tttt and Tip, events, at least one pion track is well re- 
constructed, so that its momentum in the Kg rest frame peaks around 
206 MeV, as expected for Kg — * 7i + tt~ decays. The signal distribution pop- 
ulates higher values of p* (see Fig. 5). Therefore we require for both tracks 
p* > 220MeV/c and p* n (+) + p*(-) > 478MeV/c, thus rejecting ~ 99.9% of 
tttt and 7iyx events and ~ 8% of signal events lying in the low M ee tail. 

To further reject 37r events we follow a two-step procedure. First, we exploit 
the fact that in ~ 65% of the cases the two photons originated in n decay are 
observed. Each 7 cluster is identified by TOF through the requirement 8t = 
t c i — r d/c < 5cr t . The number of detected 7's is shown in the left panel of Fig. 6 
for MC signal and 3n events. We reject events with A 7 > 2, thus introducing a 
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Fig. 5. Scatter plot of track momenta p%{+) and p%{—) (assuming it mass) in the 
Ks rest frame for MC background (left) and MC signal (right); selected region is 
shown by the solid line. 

~ 0.1% loss of signal eventsQ Second, we cut on the total missing momentum 
of the decay, evaluated as p m i SS = \p</> — Pl — Ps\, where pl,s are the neutral 
kaon momenta, and p<j> is the momentum. The Ks momentum is evaluated 
from the charged track momenta, while the Kl momentum is measured from 
the -ft'crash cluster position and the boost. The p miss distribution is shown in 
the right panel of Fig. 6 for MC signal and 3n events. We reject events with 
Pmiss > 40 MeV, thus introducing a ~ 2% loss of signal events while reducing 
the 3ir background to a negligible level. 



> 



18000 
16000 
14000 
12000 
10000 
8000 
6000 
4000 
2000 




MC 3jt 
□ MC signal 



tS 10 3 



10' 



10 



MC 3ji 
□ MC signal 




100 200 300 400 500 600 

p ■ (MeV) 

1 miss v ' 



Fig. 6. Distribution of iV 7 (left) and p m i ss (right) for MC signal (open histogram) 
and background events (gray histogram). 



A comparison between data and the expected background, N^9, in the M, 



Rejecting events with A^ 7 = 1 would improve background rejection, but would 
also introduce a systematic error in the evaluation of signal efficiency related to 
photon radiation in the final state. 



S 



sidebands after p* and iV 7 cuts is shown in Table 1. This proves the reliability 
of the background simulation and of the normalization procedure. No event is 



Cut 


Region 1 
Data ATMC A/a 


Region 3 
Data A/a 


p* > 220 MeV 

iV 7 < 2 


3738 3980(100) -2.09 
1516 1720(60) -2.83 


12107 12140(230) -0.11 
5041 5090(120) -0.36 



Table 1 

Data counts and MC background estimates in the M ee sidebands after each cut; the 
difference A between data and MC in units of the total error is also reported. 

observed in region 1, while one event remains in region 3 after p m i SS cut, for 
both data and MC background. The surviving MC event is a 37r decay. 



At this stage of the analysis, in region 2 we count three events for data and 
three hit events for MC background. To improve background rejection, we 
exploit the particle identification capability of the calorimeter. For this pur- 
pose, we evaluate the difference St(e) = t c i — Lj ' (3{m e )c between the measured 
cluster time and the expected particle time of flight under mass hypothesis 
m e , where P(m e ) = p/{m 2 e + p 2 ) 1 ^ 2 . The high rejection capability provided by 
the TOF is demonstrated in Fig. 7, where the scatter plot for St(e) of the two 
tracks is shown for MC signal and background events, before application of 
the p*, Nry and p m i SS cuts. The best rejection is obtained by using the sum 
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Fig. 7. Scatter plot of St for the two charged tracks for MC signal (box) and back- 
ground events (gray), before background rejection cuts are applied. 

S5t = 5t(e + ) + St(e~). The scatter plot of S5t as a function of M ee is shown in 
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Fig. 8 for all of the events surviving the background rejection cuts. The signal 
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Fig. 8. Scatter plot of SSt versus M ee for data (■), MC signal (gray), K$ 
(o) and 4> vr + 7r~7r (□) after background rejection cuts. 
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box is then defined as follows: 



477 < M ee < 510 MeV, 
-1.2 < SSt < 2.1ns. 



(6) 



This corresponds to a ~ 9a cut on SSt. The lower bound of the M ee range 
has been set to clearly define the fraction of IB spectrum which is accepted 
in our selection: all K$ — > e + e~ events with a radiated photon with en- 
ergy above 20 MeV are rejected. This cut reduces to a negligible level the 
other contribution to radiative decay, K$ — > 77* — > 7e + e~, which is strongly 
peaked for M ee ~ 2m e . Following Ref. [9], we evaluated BR(Ks — > 77* — > 
7e + e~, E* < 20 MeV) ~ 6 x 10 -12 , which is far beyond our experimental sen- 
sitivity. After the signal box cut we evaluate a signal efficiency given the tag 
e sig (sele) = 0.479(6). 

Applying the signal box selection, we observe no event on data. Equally, no 
MC background event falls in the signal box, see Fig. 8. In the conservative 
assumption of no background, we obtain the upper limit on the expected 
number of signal events UL(iV sig ) = 2.3, at 90% CL. 
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4 Systematic uncertainties 



Since no background subtraction has been made, there is no need to asses any 
systematic error on the scale factors applied to the MC. The selection efficiency 
for Ks —>■ e + e~ has been corrected to take account of small differences between 
data and MC tracking efficiency. The latter has been evaluated on a <p — > 
7r + 7r~7r°control sample, both for data and MC. A systematic uncertainty of 
0.9% has been evaluated by varying the correction in its allowed range. In 
order to evaluate the systematic uncertainty related to the cuts on AM ee , iV 7 
and Pmiss, we have used a K s — > ty + ti^ control sample, selected by a tight 
cut around the kaon mass for both data and MC events. The systematics 
on the previous cuts has been then evaluated on the control sample as the 
difference between data and MC efficiencies for each of the above requirements, 
as listed in Table 2. Finally, we checked the effect of IB photon emission on the 
selection efficiency. The presence of radiated photons affects indeed the shape 
of the momentum distribution, and thus the efficiency of the p* selection. 
A systematic uncertainty of 0.8% has been evaluated by comparing results 
obtained with or without the inclusion of photon radiation in the final state. 
All of the contributions to the systematic uncertainty are listed in Table 2. 



Source 


Fractional error 


Tracking 


0.£ 


% 


AM ee 


1.4 


{% 


iV 7 


0.1 


% 


Pmiss 


1.3% 




Oi 


% 


total 


2.2 


% 



Table 2 

Summary of systematic uncertainties. 
The total fractional error is 2.3%. 



5 Upper limit evaluation 

The upper limit on BR(K S — > e + e~) is evaluated by normalizing UL(A r sig ) to 
the number of Ks — > 7r + 7r~ events, N n + W -, counted in the same sample of Ks 
tagged events: 
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UL (BR(K S -> e+e~)) = 

Ctt+tt- (^crash) (sele) BR{K S -> 7T+7T") 

UM iV sigJ x 7T7 r~ X 7 — n~ X TT 

e s ig(A crash ) e sig (sele) AV+t,- 

where e(A cras h) and e(sele) are the tagging efficiency and the selection ef- 
ficiency, given -fr cr ash tag, for each channel. K$ — > tt + 7t~ events are iden- 
tified by requiring the presence of two tracks of opposite charge, selected 
with the same cuts as for Kg — > e + e~, with no additional requirements on 
invariant mass, kinematical quantities, and particle identification. The se- 
lection efficiency for both channels is evaluated from MC, with corrections 
extracted from data control samples. We obtain e n +^- (sele) = 0.6102(5) and 
e S i g (sele) = 0.479(6). The ratio of tagging efficiencies slightly differs from unity, 
e^+^(A cras h)/e sig (A cras h) = 0.9634(1). This dependence of the tagging effi- 
ciency on the Kg decay mode is due to a small difference in the determination 
of the event-To in presence of electrons or pions in the final state, which affects 
the measurement of the Kl velocity. This bias is evaluated from data using 
K s -> 7T+7T- and K s -> 7r°7r° events [10]. Using iV>^- = 217422 768 and 
BR(A's — > 7r + 7r~) from Ref. [2], we obtain: 

UL (BR (k s -> e+e" (7))) = 9 x 10~ 9 , at 90% CL. (7) 

The effect of systematic uncertainty (see Sec. 4) on the BR evaluation is 
accounted for by a Gaussian smearing of the total efficiency in the UL cal- 
culation. Our measurement improves by a factor of ~ 16 on the CPLEAR 
result [3], for the first time including radiative corrections in the evaluation of 
the upper limit. 
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